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I. INTRODUCTION 

Current interest in studies of relativistic heavy nucleus 
collisions is based on the expectation that fundamentally im- 
portant physical phenomena may occur as a result of the 
formation of high-density, high-temperature nuclear matter. 
Under such extreme conditions, matter may undergo a tran- 
sition into a deconfined quark-gluon plasma phase [1], The 
required conditions may have existed in the early universe, 
and they may be created in the interiors of neutron stars and 
in central collisions of energetic heavy ions. This last possi- 
bility provides an opportunity to study such extreme condi- 
tions in terrestrial laboratories. If high-multiplicity lead-lead 
central collisions are characterized by sufficiently high trans- 
verse momenta p, n and central pseudorapidity densities 
dN/dij, the energy densities may reach the level at which a 
quark-gluon plasma (QGP) could be formed [2], Although 
the produced particle multiplicities and their space angle dis- 
tributions will surely be dominated by common features that 
reflect kinematical constraints and variations in the impact 
parameter, new phenomena (if they exist) may be observable 
above this anticipated background in forms such as very 
large multiplicities, nonstatistical variations, or fluctuations 
in the distributions of the secondary particles. 

In December 1994, 208 Pb ions were accelerated at CERN 
to a momentum of 158 GeV/c per nucleon, by far the 
highest-energy ultraheavy nucleus beam ever produced. 

The Krakow-Louisiana-Minnesota-Moscow collaboration 
(KLMM, CERN experiment EMU- 13) exposed a series of 
nuclear emulsion chambers with Pb targets to this beam in 
order to study charged particle multiplicities and angular dis- 
tributions from interactions in the symmetric lead-lead sys- 
tem. Emulsion’s excellent spatial resolution allows accurate 
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track counting and angular measurement, with relatively 
small systematic uncertainties. In this paper, we present the 
first results from the measurement of a sample of 40 of the 
highest-multiplicity Pb-Pb collisions. In this analysis we 
consider only the gross properties of the angular distributions 
and the multiplicities. However, individual event multiplici- 
ties are sufficiently high in these collisions that it is now 
possible to search individual events for deviations from the 
behavior expected from models based on incoherent super- 
positions of nucleon-nucleon collisions. The investigation of 
fluctuations and the study of individual events will be sub- 
jects of a future study. 

II. EXPERIMENT AND ANALYSIS PROCEDURES 

The emulsions were exposed perpendicular to the beam in 
chambers of 20 emulsion plates each, spaced out over a dis- 
tance of approximately 1 7 cm from the first Pb target to the 
final emulsion plate. Each emulsion plate consisted of a 
200 fim thick acrylic base coated with a 55 ^ini Fuji ET7B 
emulsion layer on each side. An emulsion chamber is an 
extremely “light” detector, as each plate consists of only 
— 0.06 g ! cm of material. (Most tracks are measured before 
they pass through four such plates.) Each of the 32 chambers 
had a 10 cm X5 cm front area, and held three to four 
100 gm thick lead target foils. The exposure of the cham- 
bers to the beam resulted in an average of —350 primary Pb 
ions/cm across the face of the chambers, concentrated in 
three 1.5X2 cm 2 beam spots. This density was small enough 
to ensure a low delta-ray background and to keep to an ac- 
ceptably low level the number of events cut because a non- 
interacting primary was too close. 

To select a sample of relatively central interactions, the 
emulsion plates directly below each target were visually 
scanned for high-multiplicity events. After the initial scan- 
ning selections were made, each event was examined in all 
the plates upstream of the interaction and rejected if the pri- 
mary was noticeably less ionizing (approximately five charge 
units) than nearby Pb tracks or if the primary had suffered an 
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additional observable interaction. The plates adjacent to the 
target allowed rejection of interactions occumng in emulsion 
rather than in the lead target. The event was also examine 
downstream and rejected if the remnants of the project! e 
contained fragments noticeably heavier than alphas. (On y 
two events were rejected on this basis.) Events with nearby 
(60 urn) noninteracting primaries which might obscure sec- 
ondary tracks were also rejected. These high-multiplic.ty 
events are as conspicuous in the emulsion as the Pb pnma- 
ries themselves. Few if any of the very largest events are 
missed in scanning. However, the appraisal of multiplicity 
during scanning is very rough, and therefore we expect a 
gradual roll-off of scanning efficiency at low multiplicities^ 
Events with charge multiplicities above -1000 are scanned 
efficiently, but those with lower multiplicities are sampled 
incompletely. The smallest scanned event has a multiplicity 
of 590. Scanning efficiency is discussed further in Sec. in 
connection with the multiplicity distribution. 

As a result of the selection process, we have chosen 
events for analysis at a rate of ( l.42±0.18)X 10 event per 
incoming primary. By using the parametnzation of the 
charge-changing cross section for ultraheavy ion interactions 
found by Nilsen cl al. [3] and Geer cl al. [4], we expect a 
nuclear charge-changing cross section for 158 GeV/c per 
nucleon Pb-Pb interactions of 6.9 b. Using this calculate 
cross section, we estimate that we have selected 
( 22 . 2 ± 2 . 7 )% of all nuclear charge-changing interactions in 
the lead targets of the scanned chambers. 

To distinguish individual produced particle tracks emanat- 
ing from a common vertex (i.e.. the desired signal) from 
various backgrounds (delta rays, emulsion fog, emulsion sur- 
face imperfections, and particles from other events), one 
needs an image with micrometer resolution or better in all 
three dimensions, including depth. A charge-coupled device 
(CCD) camera-equipped microscope with stepper motors 
controlling all three microscope stage axes is used for this 
analysis. The acquisition is controlled by software which 
steps the focus vertically in 0.8 /rm increments through the 
emulsion layer and automatically detects the surfaces of the 
emulsion to begin and end acquisition. Depending upon the 
exact emulsion thickness, approximately 20 frames are ac- 
quired in each focus sequence. The image analysis software 
[51 searches the focus sequence for a persistent senes of dar 
pixels radiating out from a common vertex, while rejecting 
isolated dark grains and tracks which do not point back to the 
vertex The track ••darkness,” a measure of the ionization 
density is also recorded in order to distinguish minimum 
ionizing tracks from those of alphas and heavier projectile 

^Projectile fragments are expected to be confined to the 
very forward direction. Figure 1(a) relates the track darkness 
to the track emission angle 6, and shows a population o 
dark fragments mostly confined to a 2 mrad cone igure 
1 (b) shows the darkness distribution for individual tracks 
inside the forward 2 mrad cone, corresponding to pseudora- 
pidity «=-ln[tan(0/2)] = 6.9. Two peaks can be seen cor- 
responding to minimum ionizing particles and to heavier par- 
ticles (mostly alphas). We have identified tracks within this 2 
mrad cone with darkness less than 15 as minimum ionizing 
particles and tracks with darkness of 15 or more ^ frag- 
ments. The rms opening angle of the particles identified as 
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FIG. 1. (a) Track darkness vs opening angle 6. (b) Darkness 
distribution of all tracks in the forward 2 mrad cone. 

fragments is -0.7 mrad ( 77=8.0). 

Fields 108 /zmX 140 /zm across are digitized in an aver- 
age of nine plates along the axis of the event, and successive 
measurements from the individual plate sides are then fitted 
together to reconstruct the tracks in the event. By comparing 
the reconstructed tracks to their constituent measurements, 
we have determined the imaging system’s pair resolution to 
be 1 .0 am and the rms scatter of individual measurements 
within an emulsion layer to be 0.2 gm. To further discrimi- 
nate secondary tracks from backgrounds, measurements in at 
least two emulsion layers are required within - 1.0 /zm ot 
each track [5]. This requirement results in the suppression ot 
tracks below 77 = 2 . 6 . All tracks in the data sample have been 
fully measured inside the 775 * 2.9 cone. In each event, the 
track detection efficiency and background rejection are esti- 
mated for each measured emulsion layer by counting the 
missing and rejected measurements in the successive plate 
sides, respectively. The image processing software detects 
tracks with an average 96% efficiency or better for 772 * 2 . 6 . 
We have compared a sample of events reconstructed me 
software track by track with manual measurements. These 

comparisons agree to within 5%. 

In the transverse plane, the fitted track location has a sta- 
tistical uncertainty of -0.2 fxm and tracks typically leave 
the field of view at transverse distances -40 /zm from the 
event axis; the resulting 0.5% uncertainty in the transverse 
position corresponds to £ 77 = 0.005. A systematic uncertainty 
in the transverse positions derives from the absolute ^ etei ™' 
nation of the event axis. This is measured manually under the 
microscope by observing the positions of nearby noninteract- 
ing primary ions as reference tracks. The reference track po- 
sitions are determined to 5 gm; over a typical distance of 
3.3 cm (corresponding to 15 emulsion plates), this results in 
a typical systematic uncertainty of 0.15 mrad in the absolute 
positioning of the event with respect to the reference system. 
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FIG. 2. (a) Pseudorapidity distribution of the highest- 

multiplicity measured event (solid line), and the average of nine 
simulated fritiof events with similar multiplicities (dotted line) (b) 
The mean pseudorapidity distribution for the entire measured data 
sample (solid line) and that for a set of fritiof events selected with 
the same multiplicity distribution as the data (dotted line). Inset 
shows the region above 77=6 in detail. 

The uncertainty in the longitudinal track positions has a sta- 
tistical component which is greatest at large angles but does 
not exceed 1%, and an estimated 1% systematic component 
due to uncertainties in the absolute mechanical spacing be- 
tween plates during the exposure. The overall uncertainty in 
the pseudorapidity ranges from -0.01 at small 77 to 0.03 at 

17-6. The value of the pseudorapidity loses significance be- 
yond r]= 9. 


III. PSEUDORAPIDITY DISTRIBUTIONS 

Figure 2(a) shows the pseudorapidity distribution for the 
single event with the highest multiplicity. In order to com- 
pare the data to expectations based on an incoherent super- 
position model, we have simulated a sample of 1267 
Pb-Pb collisions using the fritiof 7.02 Monte Carlo code 
16J with an unrestricted range of impact parameters. In this 
preliminary study we have run fritiof in its default configu- 
ration. The dotted curve shows the average pseudorapidity 
distribution of the nine simulated events with restricted mul- 
tiplicities (V2.9-6 within the region 2.9=s tf<6 which most 
closely match that of the measured event. (We base our win- 
dow on the region above 17=2.9, where all tracks are mea- 
sured in at least two layers, and below 77=6.0, above which 
spectators are expected to appear in the measured data. Indi- 
vidual spectators are not included in the fritiof 



TABLE I. Central and semicentral data 

sets. 

Sample 

No. events 

<*z» 2) 

<*prod> 

Central 

21 

0.9±0.8 

1 3 1 4=t 210 

Semicentral 

19 

5.3±2.2 

845 ± 160 


calculations.) The two distributions are in good agreement 
(Xv 0.83 over the entire measured range). 

The mean pseudorapidity distribution (dN/dr)) for the 
entire data sample of 40 events is shown in Fig. 2(b) as the 
solid line. We have matched the measured events with 40 
events selected from the fritiof set with restricted multi- 
plicities most nearly equal to those of the real events and 
have plotted their average distribution as the dotted line. In 
t e region between 2.9 and 6, the difference between the 
distributions in Fig. 2(b) corresponds to a \ 2 per degree of 
freedom of 1.33. Again, the shapes of the distributions agree 
well except for 77>6.5, where the data show the expected 
contribution by spectators. We note that in Fig. 2(a) the mea- 
sured and calculated distributions agree well even up to the 
highest 77 values. In this most central event, few if any spec- 
tators are observed. 

To study the shape of the pseudorapidity distribution in 
the forward cone, we have separated the data set into a “cen- 
tra sample of 21 events containing two or fewer projectile 
fragments and a “semicentral” sample comprised of the re- 
maining 19 events (Table I). To compare the shapes of the 
distributions, we have normalized their areas between 
77-2.9 and 77=6 to the area of the mean inclusive distribu- 
tion, and plotted the normalized distributions for t )>6 in 
Fig. 3 along with the fritiof distribution from Fig. 2(b). The 
semicentral sample shows a component above the (spectator- 
free) fritiof prediction which is almost completely absent in 
the central sample, suggesting that fritiof predicts the shape 
of the produced forward distribution reasonably well, and 
that the “central” sample consists of events in which almost 
all of the projectile nucleons participate in the interaction. 



FIG. 3. Comparison of the shapes of the V >6 region for central, 
semicentral, and the spectator-free fritiof distributions. See text for 
details. 
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FIG 4 Relationships of forward charge and multiplicity density 
to multiplicity for several regions of the pseudorapidity distnbu- 
tions (a) Produced multiplicities in three intervals. Fits are con- 
strained to pass through the origin, (b) Total charge in several for- 
ward cones. The fits are all statistically weighted. 

(According to FRITIOF, our central event sample should con- 
tain an average of 16 spectator protons distributed over the 
pseudorapidity range 775 = 6 . We would therefore expect to 
see an excess of the number of measured tracks above the 
value for the FRITIOF pions equal to 16. In fact, we see an 
excess of 2.4± 5.0. These values differ by 2.1a, suggesting 
that FRITIOF may be overestimating the pion production in 

the forward direction by perhaps 30%.) 

Deviations from the superposition model, if they occur, 
might be expected to be strongest in the largest events. We 
look for trends in pseudorapidity shape with changing event 
size in Fig. 4. In Fig. 4(a) the mean pseudorapidity density 
appears to be directly proportional to the restricted multiplic- 
ity both near the peak and in the forward direction. In par- 
ticular, there is no indication of a flattening of the centra 
peak even for the high-multiplicity central events. This scal- 
ing implies that on average the shapes of the pseudorapi ity 
distributions are independent of the event multiplicity. This 
linear behavior is reproduced by fritiof. Table II compares 
the one-parameter linear tits shown in Fig. 4 and the corre- 
sponding fits to the FRITIOF data. The shapes of the pseudo- 
rapidity distributions agree quantitatively with FRITIOF from 
n =2 9 to 17 = 6.0 up to the highest measured multiplicities. 
(The^ uncertainties in Table II are statistical only. The 5% 
difference between the measured and the calculated slopes in 

TABLE II. Rate of increase of pseudorapidity densities with 


interval 

2.9-3.6 

4-5 

5-6 

Data 

FRITIOF 

0.457 ±0.004 
0.448 ± 0.001 

0.340± 0.004 
0.337±0.001 

0.1 58± 0.003 
0.1 67 ±0.001 


the ^= 5-6 interval is on the same order as our systematic 
counting error, and does not appear to be significant.) 

Figure 4(b) shows the total (unsigned) charge in live 
cones centered on the beam axis from i 7 > 5.5 to i?>8.0. The 
fits are shown for all five cones. For clarity, the data from 
only three representative cones are shown. The forward 
cones include spectator protons and fragments as well as 
some produced particles. We have assumed that the frag- 
ments are all alphas, and calculated the total charge in the 
interval accordingly. As Figs. 2 and 3 illustrate, the region 
forward of 77 = 6.5 contains most of the spectator contribu- 
tion. In this region, increasing the multiplicity (and the cen- 
trality) decreases the number of spectators and therefore the 
total forward charge. Widening the cone to include 77 > 5.5 
includes enough produced particles so that the charge in this 
cone increases with increasing multiplicity. 

The data in Fig. 4(b) are consistent with a linear relation- 
ship between event multiplicity and total charge. An addi- 
tional test of linearity is possible for the ? 7 > 5.5 and 
? 7>6 0 cones, which contain essentially all of the spectator 
charge. Very peripheral events must therefore have a charge 
of nearly 82 inside these cones. This is what the linear ex- 
trapolations predict. Larger cones have charge intercepts that 
are also consistent with 82. 

Summarizing, the pseudorapidity distributions are consis- 
tent with superposition in general, and agree well wit 
FRITIOF in particular. The shapes of the distributions are in- 
dependent of multiplicity. When we compare the shapes ot 
the measured Pb-Pb distributions to the shapes of simulated 
events with similar multiplicities, we see no significant di - 
ferences except those in the forward region, which can be 
attributed to spectators. 

IV. MULTIPLICITIES 

To estimate the produced charged particle multiplicities 
(i.e., the multiplicity excluding spectators) over all angles, 
we have scaled the restricted multiplicity N 2 9 - 6 b y a factor 
N nH IN 2 9-6 = L82 ± 0.06 determined from the FRITIOF 
sample with N ptod > 600 (to mimic our scanning selections). 
Adding the uncertainty in the scaling factor in quadrature 
with the estimated systematic uncertainty based on our com- 
parisons of manual and automated reconstructions, we esti- 
mate a typical uncertainty in the produced multiplicity of 
6 %. The produced multiplicity for the largest event [Fig. 

2(a)] is then 1729± 100. . 

The multiplicity distribution of our 40 measured events is 
presented in Fig. 5. We estimate in Sec. II that we have 
analyzed (22.2±2.7)% of all events in the chambers, lo 
make a direct comparison with the data, we calculate the 
FRITIOF multiplicities using the same prescription 
,y _j 82 A/ 2 9-6 as uset * to estimate the produced multi- 
plicities of the measured events. ( FRITIOF multiplicities com- 
puted using the entire 77 range produce a distribution which 
is very similar to the one shown, but which falls off some- 
what more steeply around 1850.) As expected from our event 
selection technique, we appear to undersample events with 
multiplicities less than 1000. At higher multiplicities, there is 
no evidence for an enhanced production probability. Indeed, 
we see fewer events above V prod = 1400 than expected. This 
apparent deficit is statistically unconvincing, but intriguing. 
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Produced particle multiplicity N pro<j 


FIG. 5. Probability distribution dPldN proli of the estimated pro- 
duced part.de multiplicity ^=1.82/^. The distribution of 
the data (solid line) has been normalized to an area of 0.222 based 
on the calculated cross section and event selection efficiencies The 
dotted line shows the results from an unbiased fritiof sample nor- 
malized to an area of unity. The shaded region shows the central 
events with two or fewer fragments. The right-hand axis shows the 
number of events in each multiplicity bin. 

It cannot be fully explained by normalization uncertainties in 
"prod or our scanning rate. 

To further investigate this possible deficit of large events, 
we examine the spectator region in greater detail. As the 
impact parameter b decreases to 0 , the number of spectators 
decreases. By using fritiof to estimate the number of pro- 
duced particles in the forward region, we can calculate the 
multiplicity N 0 corresponding to events with no spectators, 
i.e., events in which the forward multiplicity is entirely due 
to produced particles. This multiplicity turns out to be rather 
insensitive to the fritiof model assumptions. Figure 6 (a) 
shows the total charge Z, >6 in the cone ij> 6.0 vs N d 
The solid circles represent the central sample with two°or 
fewer fragments, and the large open circles represent the 

!rtKi 1 C f mral sampIe with more than two fragments. Figure 
4(b) shows that this cone contains essentially all of the spec- 
tator charge. The total charge of the fritiof events inside the 
? 7 > 6.0 cone has therefore been calculated by adding the pro- 
duced forward multiplicity (the “pion base line,” shown as 
small crosses) to the spectator charge, (fritiof does not 
propagate individual spectators, but does report the total 
spectator charge.) The fritiof calculation of Z >6 is dis- 
played as the small points in the top band. The fritiof dis- 
tribution converges to charge 82 on the left, and merges into 
the pion base line near " prod = 1850, which a zero impact 
parameter (b- 0 ) run confirms as the mean multiplicity N n 
oi head-on events predicted by fritiof. 

The fritiof distribution lies significantly above the mea- 
sured points. In addition, the Z r;>7 and Z, >8 distributions in 
Figs. 6 (b) and 6 (c) merge into the pion base lines near 
"prod 7 1500, not near the expected /V prod = 1850. We cannot 
explain the difference as a systematic counting error or in 
terms of a bias introduced by our event selection criteria. The 



N prod 


FIG. 6. Measured charge in several forward cones (large solid 
circles are the central sample; large open circles are the semicentra] 
sample), predicted produced charge (small crosses), and predicted 
total charge (points in upper band). Note that the vertical scales 
vary. The straight lines are statistically weighted fits, with the pion 
base line fit constrained to pass through the origin. 

intercept of the fit to the measured events at Z >6 = 82±4 
argues against a large systematic error in fragment charge 
assignments. In any case, such an error would not greatly 
effect the central sample, which has an average of only 0 9 

reaf mentS PCr CVent WC C ° nclude that the disc repancy is 

The difference indicates that fritiof cannot be correctly 
predicting both " 0 and the pion base line in the forward 
irection. We first consider the possibility that fritiof pre- 
dicts N 0 correctly and that the difference is entirely due to an 
incorrect pion base line. If " 0 = 1850 as fritiof predicts one 
consequence is that our so-called “central” sample is not 
actually very central, despite the relative lack of alphas and 
heavier fragments. From Table I, we estimate that these 
events would have on average 82X(1850- 1 3 14)/1850= 24 
spectator protons, equal to the entire mean multiplicity for- 
ward of 17 = 6.5 (25± 1 ). Thus, the produced particle pseu- 
orapidity distribution, which agrees with fritiof to within 
5% up to 17 = 6 . 0 , would have to abruptly cut off around 
r?-6.5, and the tracks forward of 17 = 6.5 in Fig. 3 would 
have to be essentially all spectators. Figures 6 (b) and 6 (c) 
confirm that for the data to be consistent with N 0 = 1850, an 
essentially complete absence of produced particles is re- 
quired in these cones. The agreement in Fig. 3, the deficit in 
the multiplicity distribution, and the lack of fragments in the 
central sample all favor the interpretation that the difference 
between the data and fritiof in Fig. 6 (a) is not entirely due 



53 


HIGH-MULTIPLICITY LEAD-LEAD INTERACTIONS AT 158 


3049 


TABLE III. Spectator depletion analysis in five forward cones. 


Cone N 0 0 ^ 


J syst 


Nk-N o 


N ’ 0 Sensitivity 


6.0 1370 

6.5 1430 

7.0 1480 

7.5 1470 

8.0 1570 


60 

80 

50 

60 

110 


70 

70 

70 

70 

80 


330 

200 

120 

80 

60 


1700 ±340 1.49 

1620± 220 0.59 

1600± 150 0.31 

1550± 120 0.17 

1 630 ± 150 0.13 


to an incorrect model of the forward region, but is in large 
measure caused by FRlTlOF’s overestimate of produced mul- 

tiP We now consider the case in which FRITIOF correctly 
models the forward pseudorapidity distribution, but overesti- 
mates the produced multiplicities. In this case, it is possible 
to estimate the number of spectators in the measured ever i s 
by subtracting the pion base line. The mean multiplicity N 0 
of head-on events, which have almost no spectators is then 
estimated by the intersection in Fig. 6 of the fits to the mea- 
sured events and to the pion base line. (For simplicity, we 
neglect the small correction due to the fact that even head-on 
events probably have an estimated four charged spectator 
nucleons. This causes us to slightly overestimate V 0 _ ) 
FRITIOF’ s total charge distribution in Fig. 6(a) crosses the 
produced particle line at 1840, which agrees well with the 
direct calculation of N 0 = 1850, demonstrating the reliability 
of the analysis technique. The analysis has been applied in 
five cones from t?>6.0 to r?>8.0, and the results are sum- 
marized in the second column of Table III. 

As discussed in Sec. Ill, FRITIOF may actually overesti- 
mate the forward production by an amount on the order ot 
30%, in which case the pion base line slopes in Fig. 6 are too 
steep, and the values of N 0 calculated in Table III are slightly 
too low Table III also gives a corrected value N 0 ot the 
head-on multiplicity in the case where the slope m of the 
pion base line is decreased by 30%. 

The systematic error <r syst in N 0 is dominated by the un- 
certainty in N prod . There is also an uncertainty in the frag- 
ment charge assignment which propagates into A„, but this 
contribution turns out to be negligible. Even assuming t a 
the fragments are all carbon only changes the value ot N 0 by 
40: The intersection is mainly determined by the fragment- 
poor central points near the pion base line. The systematic 
uncertainty due to the uncertainty in the forward production 
can be estimated from N' 0 -N 0 . Combining the statistical 
and systematic uncertainties, we find values of N 0 from the 
five cones ranging from 1550± 120 to 1700± 340, all smaller 
than the FRITIOF value of 1850. 

The sensitivity 


spending to no pions at all in the V>?-5 cone In this case, 
the head-on multiplicity would rise to 1730, still below the 
FRITIOF value of 1850. 

Finally, we combine our multiplicity and pseudorapidity 
results to find a relationship between multiplicity and peak 
pseudorapidity density. We use the data from the 

2.9— 3.6 interval in Fig. 4 along with the factor ot l.o 
from FRITIOF to quantify the relationship between produced 
multiplicity and (dNIdr,) at the peak of the pseudorapidity 
distribution The best linear fit (constrained to pass through 
the origin) gives (dN p , (vi /d t;) ^ = ( 0.25 ± 0.01 ) X . 

Thus the mean pseudorapidity density for b- 0 events 
should be 390±30. The highest pseudorapidity density we 
observe in a particular event is 425. 

Summarizing, our study of event multiplicities shows a 
significant difference from FRITIOF when the forward charges 
are compared to the event multiplicities. Our estimate of the 
mean multiplicity of head-on events N' 0 is 1550± 120, cor- 
responding to a mean peak pseudorapidity density of 
390 ± 30. No matter what the forward distribution, the best 
estimate of 1550 cannot increase to more than 1730, corre- 
sponding to (dNId 77)peak = 430. 

V. DISCUSSION 

Comparison of the data to FRITIOF shows good agreement 
in the pseudorapidity distributions at pseudorapidity densi- 
ties as high as 425. There is no evidence in the data tor 
flattening of the central pseudorapidity peak, even at pseudo- 
rapidity densities 6 times higher than in experiments at simi- 
lar energies (200 GeV/nucleon O and S on emulsion L7J). 
Such flattening might be expected if a quark-gluon plasma 
had been formed [2], It should be noted, however, that even 
for the largest event, with a central pseudorapidity density o 
425, assuming <p (7r > = 350 MeV/c (FRITIOF value) and an 
interaction distance 2ct=2 fm i/c, the energy density evalu- 
ated with the standard expression [1] from Bjorken s model. 


\N 0 IN 


0 


\m/m 


€ = 


i, dN 


4irA 2/ 'l 


of N a to the pion base line slope m can be reduced as shown 
in Table III by choosing a narrow cone. However, the statis- 
tical uncertainty increases as the cone is restricted. e 
choose the j?> 7.5 cone as the best compromise. Our best 
value of N' 0 is then 1550± 120. The smallest that the pion 
base line slope can possibly be is 0 (Am/m=-l), corre- 


(where A = 208 is the mass number), corresponds 1 to only 
1.1 GeV fm -3 . Although this energy density is significantly 
higher than in previous experiments at similar energy, it may 
still be below the point at which a quark-gluon plasma 
should be formed. 

Our determination of the mean multiplicity of head-on 
events, N' 0 = 1550± 120, is significantly lower than the value 
that FRITIOF predicts. It should be noted, however, that Ada- 
movich et al. [8] report FRITIOF simulations with a mean 
production rate of 7.68 particles per nucleon-nucleon colli- 
sion, implying N 0 = 208X7.68= 1600, in agreement with our 
measurements. The suggestion that these events are s™] ler 
than FRITIOF predictions has also been made by the EMU-U 
Collaboration based on the analysis of their first two events 


'There is a great deal of uncertainty in this number [1]. NA49 [9] 
uses a prescription which gives an energy density about twice as 
high as cited here, mainly because the formation times differ by a 
factor of 2. 
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[10]. The first results of the NA49 experiment [9] showed a 
peak negative particle multiplicity dN _ Id 77=230 for central 
events, indicating a charged particle multiplicity density of 
460. This is higher than our value but perhaps consistent 
with it. 

The value of N f 0 marks the beginning of the tail of the 
multiplicity distribution. In the superposition model, the 
width of this tail is determined by the width of the p-p mul- 
tiplicity distribution and the statistics of 208 independent 
nucleon-nucleon collisions, fritiof predicts the standard de- 
viation of the b = 0 multiplicity distribution to be 60. If 
Pb-Pb interactions are indeed simply the result of indepen- 
dent nucleon-nucleon interactions, then with better statistics 
one would expect to see a rather rapidly diminishing tail 
beyond Nq with a width of approximately 60. 

Although the method used to determine Nq requires a 
model of the pion base line, it has some noteworthy compen- 
satory advantages that distinguish it from other techniques. It 
does not rely on multiplicity cuts which could bias the result. 
(Figure 5 distinguishes the central and semicentral samples 
used in Sec. Ill, but the entire data set is used in the multi- 
plicity analysis.) It is insensitive to sampling biases, and does 
not require that the tail of the distribution be fully sampled. 
The result is almost independent of the absolute calibration 
of the forward charge measurement. And it can be performed 
with a small set of carefully measured events in which the 
tracks have been individually counted. 

In conclusion, charged particle multiplicities and pseudo- 
rapidity distributions have been measured by counting indi- 


vidual tracks in a sample of 40 high-multiplicity Pb-Pb col- 
lisions. The shapes of the pseudorapidity distributions are in 
good agreement with the results expected from calculations 
based on a superposition model of individual nucleon- 
nucleon collisions. Despite calculated energy densities twice 
those of previous experiments, we see no indication of QGP 
formation in the form of flattened pseudorapidity distribu- 
tions or enhanced multiplicities. Indeed, our best estimate of 
the mean multiplicity of zero impact parameter events is 
1550± 120, about 16% lower than predicted by fritiof. 
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